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Capture rates of compact objects were calculated by using a recent solution of the Fokker-Planck 
equation in energy-space, including two-body resonant effects. The fraction of compact objects 
(white dwarfs, neutron stars and stellar black holes) was estimated as a function of the luminosity 
of the galaxy from a new grid of evolutionary models. Stellar mass densities at the influence radius 
of central supermassive black holes were derived from brightness profiles obtained by Hubble Space 
Telescope observations. The present study indicates that the capture rates scale as oc Af^^^"''*, 
consequence of the fact that dwarf galaxies have denser central regions than luminous objects. If 
the mass distribution of supermassive black holes has a lower cutoff at ~ 1.4 x lO'' M0 (corresponding 
to the lowest observed supermassive black hole mass, located in M32), then 9 inspiral events are 
expected to be seen by LISA (7-8 corresponding to white dwarf captures and 1-2 to neutron star and 
stellar black hole captures) after one year of operation. However, if the mass distribution extends 
down to ~ 2 X 10^ M© , then the total number of expected events increases up to 579 (corresponding 
to ~ 274 stellar black hole captures, ~ 194 neutron star captures and ~ 111 white dwarf captures). 

PACS numbers: 



I. INTRODUCTION 

Dynamical evidences suggest that most, if not all el- 
liptical galaxies and bulges of spirals host supermassive 
black holes (SMBHs) with masses ranging from ^ 10^ 
M0 up to few 10^ Mq The most convincing exam- 

ple is the Milky Way, in which the orbital motion of stars 
around the radio source SgrA* indicates the presence of 
a massive compact object inside a radius less than 10^'^ 
pc, having a mass of ~ 2 — 3 x 10^ Mq Q. 

These SMBHs are generally embedded in dense stellar 
environments and, consequently, they may capture stars 
which either will be promptly swallowed through the 
horizon or will inspiral by emitting gravitational waves. 
Interactions between stars and SMBHs in dense galac- 
tic nuclei have been the focus of many investigations in 
the past years, since these events are potential sources 
for the Laser Interferometer Space Antenna (LISA). The 
planned LISA noise spectral density is minimized around 
frequencies ~ 3-30 mHz, setting the black hole mass in- 
terval of interest in the range ~ 10^ — 5 x 10^ M0. Main 
sequence stars, which dominate a given stellar popula- 
tion, will be disrupted by tidal forces if the black hole 
mass is less than '--^ 2 x 10® Mq, since the Roche limit 
will be inside their gravitational radius. Thus, only the 
capture of compact objects like white dwarfs, neutron 
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stars and stellar black holes may produce a gravitational 
signal relevant for LISA. 

Inside the influence radius r^/i of the SMBH, the or- 
bits are dominated by the gravitational potential of cen- 
tral hole but are perturbed by nearby stellar encounters. 
This means that the total mass under the form of stars 
inside rbh is less than the mass of the SMBH and that 
the radius of this dynamical sphere of influence can be 
estimated by equating the potential energy of the SMBH, 
GMbh/rbh, to the kinetic energy of the stars, fcr^^, j^. 
In the one hand, stellar-stellar interactions may produce 
important variations in energy and angular momentum 
per orbital period, leading to a prompt capture of the 
star if the resulting angular momentum is less than a 
critical value Jcrit = ^GMbh/c, for a Schwarzschild black 
hole. The phase-space volume defined by J < Jcrit is 
dubbed the "loss cone" 0, IE 01 and the infall of these 
stars in eccentric orbits onto the SMBH is set by the rate 
at which relaxation processes repopulate the loss cone 
orbits. On the other hand, stars in tightly bound orbits 
may undergo a "diffusion" in J-space with a small step 
size A J per orbital period. In this case, they lose angular 
momentum mainly by gravitational radiation, suscepti- 
ble to be detected by LISA. In the absence of the central 
SMBH, the diffusion-in is balanced by the diffusion-out 
of the "loss cone" and the net flow is zero, which is not 
the case when a secular decay of low angular momentum 
orbits occurs due to gravitational radiation. 

Due to the complexity of the energy and angular mo- 
mentum transfer, predicted capture rates found in the 
literature may vary by orders of magnitude. An early 
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Monte Carlo simulation of the stellar random walk in J- 
space was performed by 0, using parameters appropri- 
ate for the core of the dwarf elliptical M32, resulting in a 
event rate of 1.9 x 10~^ yr~^ if compact stars are assumed 
to represent 10% of the total stellar population. Similar 
simulations were performed by Q but for parameters ap- 
propriate to the galactic center. In this case, the capture 
rate of white dwarfs was found to be ~ 2 — 3 x 10~^ yr"^ 
and those of neutron stars and black holes about one 
order of magnitude lower. Analytical estimates of the 
capture rate were performed by [lOj, who have found, in 
particular for M32, a rate comparable to that obtained by 
Analytical and Monte Carlo methods were adopted 
by and from their approach a typical event rate of 
few xlO~^ yr~^ per galaxy was derived. 

As mentioned above, the capture rate is grossly given 
by the ratio between the number of stars inside the black 
hole influence radius and the two-body relaxation time at 
this position jSi i7|] . The relaxation time is essentially a 
measure of the timescale required for uncorrelated two- 
body interactions to change the original energy of a 
given star by Ai? ^ E. This timescale is also compa- 
rable to that required for the specific angular momen- 
tum to change by an amount comparable to the maxi- 
mum angular momentum for that energy, corresponding 
to circular orbits. In dynamical systems where the poten- 
tial induces precessing orbits, as long as the precession 
timescale (time required to change the argument of the 
periapse by tt) is larger than the orbital period, two-body 
interactions last longer and can be imagined as "corre- 
lated" . The consequences of coherent interactions were 
recently discussed by "^tio have concluded that the 
steady stellar flow may increase by almost one order of 
magnitude in comparison with a current of stars driven 
only by non-resonant relaxation. 

In this work the capture rate of compact objects by 
SMBHs present in the center of early-type (E-I-SO) and 
bulges of disky galaxies was estimated from the steady 
current solution of the Fokker-Planck equation in energy- 
space obtained by p^ . including resonant and non- 
resonant interactions. For a sample of E/SO galaxies, 
the stellar density at the influence radius was estimated 
from the brightness profile derived from Hubble Space 
Telescope (HST) observations, permitting the derivation 
of scaling relations either between the capture rate and 
the total luminosity of the "hot" stellar component of the 
galaxy or the mass of the central SMBH. The present 
analysis indicates that dwarf-compact galaxies, besides 
hosting SMBHs with masses in the range of interest for 
LISA, have the highest capture rates since their cores 
have stellar densities higher than giant ellipticals. Con- 
volving with the mass distribution of SMBHs, the ex- 
pected number of events to be detected by LISA after 
one year of operation was estimated. A conservative es- 
timate, including only galaxies brighter than M32, indi- 
cates that about 9 events are expected to be detected 
by LISA. If objects as faint as Mb = — 14 are supposed 
to host SMBHs, extending the mass integration down to 



'--^ 2 X 10^ Mq, the expected number of events increases 
dramatically to values around 580. This paper is orga- 
nized as follows: in Section II, the capture rate and the 
astrophysical parameters required for its calculation are 
discussed; in Section HI, the maximum redshift at which 
an inspiral capture can be detected by LISA is computed 
as a function of the black hole mass; in Section IV, using 
the previous results, estimates of the number of events ex- 
pected to be detected by LISA in one year of observations 
are given and finally, in Section V, our main conclusions 
are summarized. 



II. THE CAPTURE RATE 

As discussed in |l2j |. resonant relaxation may increase 
the inspiral rate by almost one order of magnitude rel- 
ative to rates predicted by assuming only uncorrelated 
two-body interactions. This is a consequence of the fact 
that symmetries in the grav itational potential restrict the 
evolution of the orbits |l3 |. increasing the time interval 
during which stars exert mutual torques. In the present 
work the results obtained by are used to estimate 
the capture rate in elliptical and bulges of spiral galaxies. 
Notice that the solution obtained by ^2 when applied to 
the galactic center predicts a capture rate of white dwarfs 
of about 1.2 X 10"'' yr~^, in agreement with the results 
derived from Monte Carlo simulations by |^. 

The capture rate of compact objects i?co, including 
both resonant and non-resonant relaxation is given by 
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where phh = p{fhh) is the stellar mass density at the in- 
fluence radius rth = 2GMhh/3alf^, A is the "Coulomb 
factor" expressing the ratio between the maximum and 
the minimum impact parameter of stellar interactions in- 
side the influence radius, abh is the line-of-sight velocity 
dispersion at the influence radius and fco is the fraction of 
compact objects of the stellar population present in cen- 
tral region of the galaxy. Q{x) is the dimensionless net 
steady rate at which stars flow to energies larger than x. 
The dimensionless energy is defined hy x = E/ a^f^ , where 
E is the binding energy of the star per unit of mass. Only 
stars with energies E > Eg^ can inspiral into the SMBH 
without being scattered to a wider orbit of lower energy. 
The critical energy Eg^^ can be estimated by equating 
the inspiral timescale due to gravitational radiation losses 
with the timescale required for collisions produce a varia- 
tion A J in the angular momentum of the order of J, |ll| . 
In this case one obtains for the dimensionless energy 



Xgnj = 2.4 X 10 



(2) 



where Mg is the black hole mass in units of 10^ Mq, (7200 
is the line-of-sight velocity dispersion in units of 200 km/s 
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and the stellar density at the influence radius is given in 
Mqpc-^. 

The dimensionless flow rate Q{x) is derived from the 
Fokker-Planck energy equation including sink terms due 
to resonant and non-resonant relaxation. The numerical 
solution derived by |l2j | depends on a factor which 
characterizes the efficiency of the resonant relaxation. 
Here the conservative value x = 1 is adopted. In order 
to facilitate the numerical computations, the solution de- 
rived by ^3 was fitted by a polynomial accurate enough 
for our purposes, valid in the range 1 < a; < 550, e.g., 



logQ(a;) = -0.4906 - 0.45451og(x) 



0.38261og2(a;) 
-0.21461og^(a;) (3) 



non degenerate stars 




neutron stars 



A. The fraction of compact stars 

In eq^ besides the stellar density, the fraction of com- 
pact objects fco is one of the astrophysical parameters 
required to compute the capture rate. For a given stellar 
population, the fraction of white dwarfs, neutron stars 
and black holes depends on the age and on the initial 
mass function (IMF). 

If the spectrophotometric indices of elliptical galaxies 
are interpreted in terms of single stellar populations, then 
ages ranging from few Gyr up to tens Gyr are obtained 
Hi 45j • Ehiptical galaxies are certainly not single stel- 
lar systems in spite of the fact that the bulk of their 
stars was formed in a short lapse of time. The build up 
of chemical elements requires successive stellar genera- 
tions, indicating that those systems are constituted by 
a population mix. Multi-population models point to a 
small age spread around 12 Gyr 0. A second point 
to be emphasized, consequence of these multi-population 
models, is that the IMF must be flatter than the usual 
Salpeter's law in more massive galaxies, in order to ex- 
plain the trend of the a-elements abundance ratio with 
the luminosity of the galaxy (l&ilTj . As a consequence, 
the relative number of neutron stars and stellar black 
holes in luminous galaxies are higher than in dwarf ob- 
jects. 

Here an up-graded version of the grid of models by 
[l^ was used to estimate the fraction of compact stars 
in elliptical galaxies of different masses. Models with an 
IMF of the form £,{m) oc m^+^ were developed and were 
required to reproduce the observed trends in the color- 
magnitude diagram, (U-V) versus My, as well as in the 
plane of the indices Mg2 and < Fe >. This was obtained 
by adjusting free parameters like the star formation effi- 
ciency and the exponent x of the IMF. Stars with masses 
m lower than 9.0 Mq give origin to a white dwarf of mass 
TT^wd — 0.116m -|- 0.455, while progenitors in the range 9- 
50 M0 were supposed to give origin to a neutron star of 
mass 1.4 M0. Above 50 M0 stars are supposed to pro- 
duce a black hole of 10 Mq. The upper limit of ^ 50 
M0 to form a black hole is consistent with data on X-ray 
pulsars present in binary systems like Wray 977 ^.nd 



FIG. 1: Theoretical fractions of non-degenerate stars and of 
compact objects as a function of the absolute magnitude of 
the galaxy 

Westerlund 1 T^. 

The expected fraction of non-degenerate stars (essen- 
tially main sequence stars), white dwarfs, neutron stars 
and black holes as a function of the absolute magnitude 
Mb of the galaxy is shown in figH In very massive galax- 
ies with Mb ~ —22.0, the fraction of white dwarfs may 
reach values as high as 16%, whereas the fraction of neu- 
tron stars and black holes attain values of about 3.2% 
and 0.4% respectively. These results will be used in our 
capture rate estimates. 

B. The stellar density 

Another important astrophysical quantity required to 
compute the inspiral rate from eq^ is the stellar den- 
sity at the influence radius. Near the center, the surface 
brightness along a given axis may be written as 

/(r) =/o + /ir2+0(/) (4) 

where r is the projected distance from the center. Ac- 
cording to jl^l , a galaxy satisfying the above relation can 
be said to have an "analytic core" and its core radius 
Rc is defined by the condition I{Rc) — |-f(0). However, 
HST data have shown that real cores are not analytic but 
display shallow power-law cusps into the resolution limit, 
which brightness can be fitted by the so-called Nuker law 
e.g., 

/(r) = 2('^-^)/"/fc(-)^[l + (Zl)"](7-/3)/a (5) 
r Tb 

where 7^, rp, a, /3 and 7 are fitting parameters. Ac- 
cording to |2Q|, essentially two types of galaxies can be 
distinguished on the basis of the aforementioned parame- 
ters: a) "core" galaxies, which have an inner logarithmic 
slope I dlogJ/dlogr |< 0.3 and b) "power law" galaxies, 
which show a fairly steep brightness profile for r < rf,. 
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The former class includes luminous objects whereas the 
later includes fainter galaxies. These properties are fun- 
damental to understand the stellar environment in which 
the central black hole is embedded. 

Once the intensity profile is known, the luminosity den- 
sity PL can be derived by inverting an Abel integral |22| . 
e-g-, 



1 f" 

Jr 



{dl{r)/dr) 



dr 



(6) 



and the stellar mass density can be obtained by multi- 
plying the luminosity density by the mass-to-luminosity 
ratio Ty — M/Ly in solar units. From our models, 



logTy = -1.148 - 0.0956My 



(7) 



The stellar density at the influence radius was com- 
puted numerically using eqs|5lEland[3 The parameters 
defining the brightness profile as well as absolute mag- 
nitudes were taken from reference (24| . corrected for a 
Hubble constant Hq — 70 kms~^Mpc~^, adopted in the 
present work. Black hole masses required to compute the 
influence radius, when available were taken from p^l2^ 
or computed from the relation 



(8) 



which is slightly different from the original fits derived 
by or 27], since here only E and SO galaxies were 
included (bulges of spirals were excluded from the fit). 

Fig[21 shows the stellar density at the influence radius 
as a function of the absolute magnitude of the galaxy. 
This plot just confirms the known fact that dwarf el- 
lipticals have central densities considerably higher than 
bright E-galaxies, which is essentially a consequence of 
"power-law" profiles present in the former objects and 
the existence of a fundamental plane, correlating struc- 
tural parameters like the central velocity dispersion, the 
effective brightness and the effective radius. As a corol- 
lary, SMHBs with masses around ~ 10^^^ Mq are al- 
ways embedded in denser stellar environments, thus hav- 
ing higher capture rates, than those surrounding SMBHs 
with masses around 10^~® Mq, which have capture rates 
one or two orders of magnitude lower. 



C. The estimated inspiral rates 

For galaxies included in the sample studied by [2^ . 
capture rates were calculated from cq^ using the frac- 
tion of compact objects, stellar density and influence 
radius estimated according to the procedure described 
above. The resulting capture rate of white dwarfs scales 
with the luminosity of the galaxy as 



.31 X 10^ 



{Lb/Lb,q) 



1.31 



(9) 
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FIG. 2: Stellar densities at the influence radius as a func- 
tion of the absolute B-magnitude for galaxies included in the 
sample by Faber et al. |24i. Absolute magnitudes were taken 
from LED A database 1251. 



whereas the scaling with the SMBH mass is 
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The capture rate of neutron stars and stellar black holes 
can be obtained by multiplying the equations, above 
respectively by 0.0511 and 0.00916. In spite of these small 
rates, the gravitational signal resulting from the capture 
either of a neutron star or a stellar black hole can be 
detected much further away than that resulting from the 
capture of a white dwarf, which makes their detection 
rate comparable, as we will demonstrate later. Notice 
that due to the square dependence on the stellar mass 
density in eq^ the capture rate decreases with increasing 
SMBH masses, since more massive holes live in less dense 
stellar environments. 

If one assumes that the relaxation timescale is com- 
parable to the mean population age, e.g., 12 Gyr, then 
the stellar density and the velocity dispersion at rb/i must 
satisfy 



0.154(72^^05 



Mopc- 



(11) 



where aio is given in km/s. On the other hand, densities 
at the influence radius derived in the previous section and 
the line-of-sight velocity dispersion are approximately re- 
lated as 



p w 1.44 X lO^^cr-^-^^'' Mqpc" 



(12) 



From these two relations, in order that < 12 Gyr, 
one must have aio < 100 km/s and/or pbh > 9.0 x 10^ 
Mqpc^"^. This means that only galaxies fainter than 
Mb = — 18 have relaxation timescales less or equal 12 
Gyr and can satisfy the steady conditions required by 
eq^ Galaxies brighter than this limit are not able to 
refill their "loss cones" within a Hubble time and prob- 
ably have a non steady inspiral flow. Therefore, we 
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will restrict the application of eq. |H1 or eq. ^| to ob- 
jects fainter than the aforementioned limit, which corre- 
sponds to galaxies hosting central black holes with masses 
around - 8 x 10*^ Mq. 

Eqs. m El are valid, in principle, only for early-type 
galaxies. However, here they will be also applied to 
bulges. Bulges constitute, in general, an extension of 
low luminosity ellipticals. They follow quite well the 
Mg2 — (JiD diagram of ellipticals [23, the SMBHs present 
in the center of bulges follow approximately the same 
— cTiD relation obeyed by early-type galaxies (but 
having probably a slightly steeper slope although more 
data are needed to confirm such a trend [2l|) and stel- 
lar densities derived at the influence radius, at least for 
two cases, are in agreement with the results displayed in 
figUl These are the bulges of the Milky Way, if we use 
the results by 30] and that of NGC 1316, if we use the 
procedure described above and data from j2j| . 



III. THE REDSHIFT SPACE PROBED BY LISA 

Estimates of the number of events expected to be de- 
tected by LISA demand a previous evaluation of the vol- 
ume of the universe probed by the detector or, in other 
words, an evaluation of the maximum redshift at which 
a given inspiral gravitational signal can be seen. 

For the inspiral up to the last stable orbit (LSO), 
theoretical templates of the waveform h{t) will proba- 
bly be available and the method of matched filtering can 
be used. In this case, the resulting signal-to-noise ratio 
(SNR) is 113 



(-)' 



Sn{v) 



where Sn{v) is the spectral density of the strain noise of 
the detector and the Fourier transform of the two polar- 
ization amplitudes, /i+ and /ix, depends on the source 
orientation, on the detector beam pattern functions -F+ 
and Fy, , on the polarization and on the source sky lo- 
cation. Averaging over source positions and orientations 
one obtains l32l: 



The spectral density of the gravitational signal is given 
by j^] 



1 G (1 + ^)2 1 dE 



2 7r2c3 



(17) 



where v and v,, are respectively the gravitational wave 
frequencies at the observer's frame and at the source re- 
lated by i^e — + z)v, (ii = (1 -I- z)r(z) is the distance 
luminosity with r[z) being the proper distance defined 

by 



J^o Jo 



dz' 
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with 



Bin,, z) = ^n^ + nrn{i + z)3 



(18) 



(19) 



where fl^ — 0.7 and — 0.3 are respectively the present 
values of the density parameters due to vacuum and mat- 
ter (baryonic -I- non-baryonic) adopted in the present 
work, in agreement with the "concordance" model (|S3l)- 
In eqEI the spectral energy density of the source in the 
quadrupole approximation, which is adequate to estimate 
the SNR obtained from matched filtering 



dE 



{Gt:) 



2/3 



2/3j^-l/3 



(20) 

and the total mass 



where the reduced mass ^ m r 
M « Mbh. 

Combining eas ll7l and ITHl into 1141 one obtains 
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where 
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(21) 



(22) 



The upper limit of the above integral corresponds to the 
frequency associated to the LSO, corrected for the ob- 
server's frame, e.g.. 
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where the response function of the detector averaged over 
sky directions and polarizations can be expressed as : 



/2T67rGAf6h 



0.004384 



Hz 



(23) 



For the lower limit we take the frequency corresponding 
to about one year before arrival at the LSO, namely. 



/f^ E F^*{^,^)F\^.^) (15) 



On the other hand, the planned LISA sensitivity is usu- 
ally given in terms of the quantity 
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where rrico is in solar units and A/g is in units of 10^ Af©. 

Following the convention adopted in the LISA com- 
munity, the detectability threshold was settled to be 
S/N = 5 and from the numerical solution of eq[2] one 
obtains, for a given compact object of mass rrico, the 
maximum redshift as function of the SMBH mass Afb^,. 



6 



the local mass distribution function 
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FIG. 3: Maximum redshift (zmax) at which a gravitational 
signal resulting from an inspiral capture can be detected by 
LISA operating as a Michelson interferometer. A signal-to- 
noise ratio equal to 5 was adopted and the WD- WD galactic 
binary confusion noise was included. 



The effective noise budget sensitivity (see ea ll6|) used in 
the calculations refers to a standard Michelson config- 
uration including the contribution of the galactic 
binary WD- WD confusion noise [s^ |3|| . The resulting 
curves for white dwarfs (m^.d = 0.7 Mq), neutron stars 
(m„s = 1-4 M0) and stellar black holes {nibh — 10 Mq) 
are shown in fig|3l The farthest inspiral signal result- 
ing from the capture of a white dwarf corresponds to a 
redshift z - 0.76 and a SMHB mass of about 3.9 x 10^ 
Mq, that resulting from the capture of a neutron star 
corresponds to 2: - 1.19 and a SMBH of 3.1 x 10^ Mq 
whereas the signal resulting from the capture of a stellar 
black hole by a SMBH of mass 1.3 x 10^ Mq can be seen 
up to z ~ 3.67. 



IV. THE EXPECTED NUMBER OF EVENTS 

Once the maximum redshift z^ax at which a given in- 
spiral signal can be detected is calculated, the expected 
number of events Tco in a time interval T can be esti- 
mated from the equation 



M2 



Ml 



Rco{Mbh)dJ\f{Mbh,z) 
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dM, 
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dz 



dMbh J a 

(25) 

In this equation, Rco{Mbh) is the capture rate of a given 
compact object as a function of the black hole mass, as 
calculated in Section II. B and the term (1+z) in the de- 
nominator corrects for time dilation as the universe ex- 
pands. Evolutionary effects which may affect the mass 
of the central black hole and/or the nearby stellar dis- 
tribution were neglected in the present approach but are 
under investigation. For the SMBH mass distribution, 
the results by j33| were used. These authors derived for 



dM, 



bh 



Mo Ml 



^ Mq Y-^^-(M,^/MoY 



bh 



(l + zf (26) 



where the parameters 6**, Mq and 7 where calculated 
for the different morphological types, using the lumi- 
nosity distribution functions by [33|. The parameter 
e = 0.8 ±0.1 depends on the slopes of the adopted 
Mbh — <yiD and L — any relations. Since the luminos- 
ity function concerns the "total" galaxy luminosity, in 
the case of bulges of spirals the derivation of the black 
hole mass distribution depends on the relative bulge- 
to-total luminosity ratio for the different morphological 
types. Since average values for the various morphological 
types have a considerable dispersion, this introduces an 
important source of uncertainty in the parameters corre- 
sponding to disky galaxies. In above equation, the con- 
servation of the number of galaxies within the comoving 
volume was included but no other evolutionary effect. 
The comoving volume in eq. 1251 is given by 



(27) 



In the second integral of eq|5S| the upper limit M2 
corresponds to the critical luminosity beyond which the 
relaxation time is larger than 12 Gyr and, consequently, 
a steady stellar current can not be established. The con- 
version of luminosities into SMBH masses was performed 
through the relation 



logMbh/MQ = -0.835 - 0.439 A/b 



(28) 



If, as we shall see, our results are practically indepen- 
dent on the assumed upper limit M2, the situation is 
quite different concerning the lower limit Mi, since the 
capture rate increases for these less massive black holes. 
The lower mass cutoff is uncertain since up today obser- 
vational searches for intermediate-mass black holes have 
been discouraging. Presently, the smallest known central 
black hole is located in the dwarf elliptical M32, having 
a mass ^ 1.4 x 10^ Mq (j2Q]). A conservative estimate of 
the number of events can be obtained if this value is taken 
as the minimum mass. From an observational point of 
view, the existence of SMBHs with masses lower than the 
above limit is still controversial. Analyses of dynamical 
data on the dwarf galaxies M33 [H ll^ and NGC 205 
PH lead to robust upper limits respectively of ~ 3 x 10'^ 
Mq and ~ 4 x 10^ Mq to the putative black holes living 
in the center of these objects. Moreover, searches for in- 
termediate mass black holes in dwarf-spheroidal satellites 
of the Milky Way have also been negative ^J. However, 
these small galaxies were captured by our Galaxy and, in 
the process, collective effects may transfer energy to their 
central black holes enough to move them far away from 
the core, escaping from detection. If direct dynamical 
evidences of SMBHs in such a mass range are missing, 
indirect signals of their presence in narrow- line Seyfert-1 
nuclei exist. Using the line width- luminosity-mass scaling 
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relation established for broad-line AGNs, SMBH masses 
were estimated by for a sample of 19 galaxies, all 
in the range 8 x 10"* up to 8 x 10^ M©. Therefore, an 
"optimistic" estimate was also performed by assuming 
that central black black holes are present at least in com- 
pact dwarf galaxies and bulges brighter than Mb = —14, 
since their cores follow quite well the A/fc/i-luminosity re- 
lation fSsI . Using eq. 1281 this luminosity corresponds to a 
SMBH of about 2 x 10^ Mq . Such a lower limit were used 
for "optimistic" estimates of the capture rate of neutron 
stars and stellar black holes, since for white dwarfs the 
adopted lower limit was slightly higher ^ 6 x 10^ M©, 
corresponding to the mass for which the tidal radius is 
comparable to the gravitational radius. The results of 
our computations indicate that in the conservative case, 
9 inspiral events are expected to be detected by LISA. In 
this case, the estimated probability to have a white dwarf 
capture is 84.5% and those for the capture of a neutron 
star or a stellar black hole are 8.0% and 7.5% respectively. 
Early-type galaxies are expected to contribute to about 
53.8% of the total number of events, Sa-|-Sb bulges con- 
tribute to about 26.9% and Sc bulges contribute to about 
19.3%. 

If the lower mass integration limit is reduced down 
to 2 X 10^ Mq, there is a dramatic increase in the ex- 
pected number of events. This occurs because the cap- 
ture rate increases for lower SMBH masses as well as the 
volume of the universe probed by LISA (see figEI)- In 
this "optimistic" case, the expected number of events is 
579. The contribution of the different compact objects 
is now rather different: stellar black holes are expected 
to represent about 47.4% of the total number of events, 
neutron stars will contribute to about 33.5% and white 
dwarfs, to about 19.1%. The contribution of the different 
morphological types will be practically the same as be- 
fore, e.g., E-I-SO will contribute to about 53.2%, Sa-|-Sb, 
to about 20.5% and Sc to about 26.3%. 



V. CONCLUSIONS 

Diffusion in energy-space of stars under the gravita- 
tional influence of a supermassive black hole was recently 
considered in [l^ . including phenomenologically a term 
in the Fokker-Planck equation, corresponding to two- 
body resonant relaxation. These results were used to 
estimate the capture rate of compact objects in a sample 
of early- type galaxies. 

From a new grid of evolutionary models for early-type 
galaxies, including a variable initial mass function, the 
fraction of white dwarfs, neutron stars and stellar black 
holes was derived as a function of the luminosity (mass) 
of the galaxy, as well as the (stellar) mass-to-luminosity 
ratio, quantities necessary to compute capture rates. The 
later was used to estimate stellar mass densities by in- 
version of brightness profiles of a sample of early-type 
objects, obtained from HST observations. The analysis 
of these results indicates that the capture rate decreases 



for high luminosity galaxies, which harbor the most mas- 
sive black holes. This is a consequence that luminous E- 
galaxies have cored density profiles whereas fainter ob- 
jects have in general "power-law" profiles, providing a 
denser stellar environment to their hosts. 

Interpretation of colors and metallicity indices of early- 
type galaxies by evolutionary models indicates mean pop- 
ulation ages of about 12 Gyr with a dispersion of about 
2 Gyr. If dynamical ages are comparable, then only sys- 
tems fainter than Mb = — 18 are able to refill their loss 
cones and maintain a steady stellar current. Such a lumi- 
nosity corresponds to a SMBH mass of ^ 8 x 10^ Mq and 
constitutes a limit for the application of steady solutions 
of the Fokker-Planck equation. 

Extending our results to bulges of spirals, face their 
similarity to E-galaxies, the expected number of events 
to be detected by LISA in one year of operation was es- 
timated. Two situations were considered: in the first, 
a conservative estimate was performed by restricting the 
integration over masses in the range defined by the afore- 
mentioned upper limit and the lower limit defined by the 
lowest SMBH mass detected up today, e.g., the one lo- 
cated in M32. A second and a more optimistic estimate, 
extends the lower mass limit down to 2 x 10^ Mq corre- 
sponding to a luminosity of Mb — —14.0. The core of 
these compact dwarf galaxies, able to host those unde- 
tected up to now black holes, follow the M — a relation 
|43l | defined by brighter galaxies and bulges, which is en- 
couraging. In this case, the resulting expected number 
of events is 579. The inclusion of SMBHs with masses 
lower than few 10^ Mq increases not only the capture 
rate but also increases the volume "seen" by the detector 
(see figlSll- In the conservative case, between 7-8 events, 
probably associated to the capture of white dwarfs and 
1-2 events, probably associated to the capture of neutron 
stars and stellar black holes are expected to be seen by 
LISA after one year of operation. Reducing the lower 
mass limit, events associated to the capture of stellar 
black holes will be dominant (~ 274), followed by those 
associated to the capture of neutron stars (~ 194) and 
white dwarfs (~ 111). These aspects associated with the 
number of events and the distribution of the different 
captured compact objects reveal the potential of gravi- 
tational waves as tool for astrophysics, in particular for 
scenarios of SMBH formation. 

However, it should be emphasized that the present esti- 
mates neglect possible evolutionary effects. If mean pop- 
ulation ages of early- type galaxies and bulges are ^12 
Gyr, this means that they were already assembled around 
z ~ 3.6, which corresponds practically to the maximum 
redshift that an inspiral gravitational signal produced by 
the capture of a stellar black hole by a SMBH of mass 

- 1.3 X 10^ M0 can be detected by LISA with a SNR 

— 5. No important variations are expected to occur in 
the central region of the galaxy as a consequence of the 
capture of compact objects in such a timescale. However, 
non-degenerate stars will diffuse-in the influence radius 
and will be destroyed by tidal forces, representing a mass 
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consumption in a timescale of 12 Gyr which may be com- 
parable to the that of a 2 x 10^ Mq black hole. This 
process will modify slowly the central density profile and 
the SMBH mass. Moreover, hosting dark halos are con- 
tinuously accreting mass |4J], affecting the mass history 



of central black hole and the structure of the galaxy as 
well. These problems are presently under investigation 
through cosmological numerical simulations and will be 
reported in a future paper. 
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